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Abstract. This paper deals with the homogenization problem for a one-dimensional parabolic PDE with random 
stationary mixing coefficients in the presence of a large zero order term. We show that under a proper choice of the 
scaling factor for the said zero order terms, the family of solutions of the studied problem converges in law, and 
describe the limit process. It should be noted that the limit dynamics remain random. 

Resume. Cet article traite de l'homogeneisation d'une equation aux derivees partielles en dimension un d'espace, 
avec des coefficients aleatoires stationnaires et melangeants, en presence d'u terme d'ordre zero fortement oscillant. 
Nous montrons qu'avec un choix convenable du facteur d'echelle de ce terme d'ordre zero, les solutions du probleme 
etudie convergent en loi, et nous decrivons le processus limite. On peut noter que la dynamique limite est elle aussi 
aleatoire. 

MSC: 74Q10 
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1. Introduction 

Our goal is to study the limit, as e — > 0, of a linear parabolic PDE of the form 



^ s) = 2toK^te) ( ^ B) + Vi c W^ (t,!,f) ' " °' x€ ; 

u e (0,x) =g(x), ieM, 

where a and c are stationary random fields, and c is centered. 
Let us recall (see [1]) that in the periodic case the equation 

d 19// x\ d \ 1 / x 



(1.1) 



dt U= 2Tx\ a \-e d-x U + ~e C ' 



This is an electronic reprint of the original article published by the Institute of Mathematical Statistics in 
Annales de l'lnstitut Henri Poincare - Probabilites et Statistiques, 2008, Vol. 44, No. 3, 519-543. This reprint 
differs from the original in pagination and typographic detail. 
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admits homogcnization under the natural condition (c) = ((•) stands for the mean value) and that the 
homogenized operator takes the form 

d 1. d 2 

m u= 2 a d^ u + cu > 

with constant a and c. 

In contrast with symmetric divergence form parabolic problems, in the presence of the lower order terms 
the asymptotic behaviour of operators with random coefficients might differ a lot from that of periodic 
operators. 

Homogcnization problem for parabolic operators whose coefficients are periodic in spatial variables and 
random stationary in time, were studied in [4, 5, 14]. It was shown that, under natural mixing assumptions 
on the coefficients, the critical rate of the potential growth is of order 1/e. In this case the limit equation is 
a stochastic PDE. 

If the oscillating potential is random stationary (statistically homogeneous) in spatial variables then the 
range of the oscillations (the power of e -1 in front of potential c) should depend on the spatial dimension. 

In this work we deal with a one-dimensional spatial variable and show that the range of oscillation should 
be of order -^=. This means that for larger powers of i the family of solutions is not tight as e — » 0, while 

for smaller powers of - the contribution of the potential is asymptotically negligible. 

It turns out that the Dirichlct forms technique which is usually quite efficient in homogenization problems, 
does not apply to problem (1.1) because one cannot prove any lower bound for the quadratic form corre- 
sponding to the operator (1.1). This is due to the fact that the problem is stated on the whole line R, and 
not on a compact interval, and that the coefficients of the operator are a.s. unbounded, see the discussion in 
Section 6. Instead we use the direct approach combining the Fcynman-Kac formula with several correctors, 
Ito calculus and martingale convergence arguments. 

The main result of the paper (see Theorem 2.2) states that under proper mixing conditions the solution 
u £ of eq. (1.1) converges in law to a random field 



u(x,t)=E(g(x + V5B t )exp^| J L y t ~ x W{dy) 



where B. and W. are independent Brownian motions, E and L\ are respectively the expectation and the 
local time related to \fd,B. , and a, c are constants. 

The interpretation of this expression is given in the last section of the paper. It is shown that the effective 
equation is not a standard SPDE but rather a parabolic PDE with random coefficients. 

Let us give an intuitive explanation of our result. The Feynman-Kac formula for the solution of eq. (1.1) 
yields 



u e (t,x) = E 



where E means expectation with respect to the law of the diffusion Xf' x , the random field c being frozen (or 
"quenched"). Under the assumptions which we shall make below, one can apply a version of the functional 
central limit theorem, which tells us that 

W £ (x) = -^ f c( V -\dy 



converges weakly towards cW(x), where W is a standard Wiener process. Now the exponent in the above 
Feynman-Kac formula reads 

Jo Jr 
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where L\' e denotes the local time at time t and point x of the diffusion process {X £ ' x }. One might expect 
that the last integral converges towards the integral of the limiting local time, with respect to the limiting 
Wiener process. This is one of the results which will be established in this paper. 

Our paper is organized as follows. In Section 2, we formulate our assumptions, and the results. In Section 3, 
we prove some weak convergence results. Section 4 is devoted to the proof of the pointwise convergence of 
the sequence u E (t,x), while Section 5 is concerned with convergence in the space of continuous functions. 
Finally in Section 6 we discuss the limiting PDE. 



2. Set up and statement of the main result 



We make the following assumptions: 

(A.l) The initial condition g belongs to L 2 (R)nC b (R). 

(A. 2) The coefficients {a(x),x G K} and {c(x),x G R} are stationary random fields defined on a probability 
space (S,A,P), and we assume that 

0<c<a(x)<C, i£l,Pa.s, (2.1) 

£c(0) = 0, ||c(0)|Uoo (J3) < oo, (2.2) 

|^c(0)c(x)|dx<oo, (2.3) 

where E denotes expectation with respect to the probability measure P. 
(A.3) Let 

T x :=<r{a(y),c(y);y<x}, T x := <r{a(y),c(y);y > x}. 

We assume that the random fields a and c are 0-mixing in the following sense. Define, for h > 0, <p(h) the 
mixing coefficient with respect to the c-algebras from above, as 

<f)(h):= sup \P(B\A) — P{B)\. 

We suppose that 

4> 1/2 {h)dh<oo. (2.4) 
Consider now the family of Dirichlet forms {E e ' a , e > 0, a G £} on L 2 (R) defined by 



£ E - a (u,v) = -J ai^-,aju'{x)v'(x)dx, 

with domain 7J 1 (R). For each e > 0,er G S there exists a unique self-adjoint operator L 6 ' 17 with domain 
X>(i7' <T ), such that 

for u£V(L E - a ),v£ L 2 (R). 

For each initial point i£R, a G £ and e > 0, there exists a continuous Markov process {Xt' cr ' x , t > 0} 
defined on some probability space (f},T,f e ' x ' a ) whose generator is L ecr and which starts at time t = from 
x. The probability may depend on the three parameters e, x, a. x will be fixed throughout this paper, 
so we drop it from now on. Note that the process {X^' x ,t > 0} is in fact defined on the probability space 
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(S x f2, A (£) J- , Q £ ) (as such, it is not a Markov process), where the probability Q £ on the product space 
S x Q is defined as 



Q £ (A) = J J P(da)P e ^(d^). 



The Feynman-Kac formula allows us to write down an explicit formula for the solution of eq. (1.1): 



u e (t, x) 



(2.5) 



where E e,cr denotes expectation with respect to P e ' <7 . 

Considering assumption (A. 2), we define the finite quantities 



E[c(0)c(x)]dx, a = 



E 



1 



a(x) 



(2.6) 



In view of Theorem 5.1 and Lemma 5.1 from [11], we may state the following theorem. 
Theorem 2.1. We have the following convergence, P a.s.: 
Xf' x ^>Xf :=x + X., 

in C([0,oo)), where X is a one- dimensional Brownian motion defined on (f2,J-,P) such that E(X t 2 ) = at, 
for any t>0. 

The main result of this paper is the following theorem. 
Theorem 2.2. Let 



u(t, x) := E 



g(x + X t )exp[ ~ / L y t - x W(dy) 



where W denotes a one- dimensional standard Brownian motion defined on the probability space (£, A, P) 
and L\ is the local time at time t and point y of the process {X t ,t > 0} defined on (i?,.F,P). 
Then u s u in law in C(K+ xl), as e — > 0. 

We introduce the notation 



V€,X 

ci ^2- ]d«. 

e Jo \ e 



The first step in the proof of Theorem 2.2 is to establish the weak convergence of the pair (X t ' ,Y t ' ), 
which is done in the next section. 



3. Weak convergence 

The main result of this section is the following theorem. 



Theorem 3.1. For each t>0, 
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weakly, as s — > 0, with 
a Jm 

where, as above, L\ is the local time at point y and time t of the Brownian motion {X t ,t > 0} defined on 
(J?,.7 r ,P), and {W y ,y€ P} is a Wiener process defined on (S,A,P), so that (X,L) and W are independent. 

Theorem 3.1 will follow easily from Propositions 3.7 and 3.10, as we shall see at the end of this section. 
Note that all we shall need in the next section is both Propositions 3.7 and 3.10. 
Let us first state a consequence of Aronson's estimate, see Lemma II. 1.2 in [16]: 

Lemma 3.2. There exists k > 0, which depends only on c and C in (2.1), such that for all e > 0, r > 0, 
pf sup \X*' x -x\ >r) <kcx P ( 



nt 



We next prove the easiest part of the above result, i.e. we give a proof of Theorem 2.1, since we shall 
need some of its details later. 

Proof of Theorem 2.1. Let {\(x),x e R} be the zero mean random process given by the formula 
Jo a W 

We note that from Birkhof's crgodic theorem (see e.g. Theorem 24.1 in [3]), 

^-^0, P a.s., as \x\ -> oo. (3.1) 

Moreover, this random process satisfies the two relations: 

(a(l + X '))'(x)=0, x £ M, (3.2) 

and 

a(x)(l + x'(x))=a, 
We now define 



z; = x?*+ex 

It follows from the Ito-Fukushima decomposition (see [7] or Theorem 0.10 in [11]) and (3.2) that (here and 
further below, M x denotes the martingale part of the process X £ - x ) 



%t — Zq 



t 

i + x 



II 



AMf'\ t>0, 



e 

hence {Z E } is P a.s. a P-martingale. Moreover, its quadratic variation is given by 
ds 



/o a{Xr/e) 
It will be proved below in Lemma 3.9 that 



{Z e ) t -~at 
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in Q e probability. It now follows from well-known results that P a.s., 
Z e =S> x + Va~B, 

where {B t ,t > 0} is a standard Brownian motion defined on the probability space {fl,T , P). 
Moreover, for all T > 0, 

sup \X^' X -Z^\^0, Q £ a.s., 

0<t<T 

consequently P a.s., 
X £ ' x ^x + Va~B 

in P law. □ 
Let # denote the solution of the ordinary differential equation: 
(a$ , )\x) = c(x), 

which is defined as follows: 

1 f x 
$'{x) = - n - / c{y)dy, 

[ . (3-3) 



o \a(z) Jo 



We let 



cW e (aO :=Vial ~W(~ 



1 r ^ 



V? Jo V £ 



dy (3.4) 



and 



F B (x) :=£ 3 / 2 <Z>( - 



c(— \ dydz 



e J a(z/e)J \e 
I 



W e (z)dz. (3.5) 



/o a(^/e) 

We first prove the following proposition. 

Proposition 3.3. The sequence of random processes {W e } converges weakly, as £ — ► 0, in i/ie space C(M), 
to a standard Wiener process {W} defined on (S,A,P). 

Proof. Denote, for x > 0, W^x) = W e (x) and Wg (a;) = W e (— x). According to assumptions (A.l), (A. 2), 
(A. 3) and the functional central limit theorem (see e.g. [2], pages 178, 179), it follows that 

(w^w!)^(w\w 2 ), 
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where {W (x),x > 0} and {W 2 (x),x > 0} are mutually independent standard Brownian motions. Finally 
we denote by {W(x),x £ R} the process defined by 

W{x) := W 1 {x), fora;>0, W{x) := W 2 (-x), for x < 0. □ 

It remains to show why Theorem 3.1 follows from Theorem 2.1 and Proposition 3.3. 
First we define, for x £ R, 

f x 1 x 
k £ (x) = , M dz, fc(x) = -. 
jo a(z/e) a 

We have the following lemma. 
Lemma 3.4. As e — » 0. 

fc e -> & m C(R),P a.s. 

Proof. Since a(-) is bounded away from zero, a _1 (-) is bounded. Hence the collection of random functions 
k £ is tight in C(R). It then suffices to show that the finite dimensional marginals converge in law to those 
of the deterministic function k. But from BirkhofFs ergodic theorem, for any x%,x%, . . . ,x n £ R, 

(k E (x 1 ),k e (x 2 ), ■ ■ .,k E (x n )) 

in P a.s., as £ — > 0. □ 

Denote by C+(R) the space of continuous and increasing functions on R, and by S the map from C(R) x 
C+(R) into C(R) defined by 



S((f,h))(x):= / f(z)dh(z). 
Jo 

We have the following lemma. 

Lemma 3.5. The mapping S is continuous, from E = C(R) x C+(R), equipped with the product of the 
locally uniform topology o/C(R) x C(R) ; into C(R), equipped with the locally uniform topology. 

Proof. It suffices to show that for each N > 0, if {(/„, h n )} C C([-N,N\) x C+([-N,N}), and 
sup (\fn{x) - f(x)\ + \h n (x) - h(x)\) -»0, asruoo, 

|x|<JV 

then 

sup \S(f„,h„)(x) - S(f,h)(x)\ ^0, asruoo. 

|i|<Af 

But this follows from Lemma 5.8 in [8]. □ 

We now have the following lemma. 
Lemma 3.6. As e — » 0, 

(W e ,F e ) — > {W,F) 

in C(R) x C(R) m Zaw, where F(x) = § J"^ dz and W e and F e are defined in (3.4) and (3.5) respectively. 
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Proof. Since k e converges to a deterministic limit, it follows from a well-known theorem (see e.g. Theo- 
rem 4.4 in [2]) that the pair (W £ , k £ ) converges. Hence from Lemma 3.5, the pair (W e ,F e ) = (W s , S(W £ , k e )) 
converges. □ 

The next step is to show that the triple (X e ,W e , F e ) converges. This is essentially a consequence of the 
three following facts: X £ converges, (W £ ,F £ ) converges, and the two limits X and (W,F) are defined on 
(fi, T, P) and (E,A,P) respectively. We now prove that fact rigorously. 

Proposition 3.7. For any as e — > 0, 

(X £ ' X ,W £ ,F £ )^(X X ,W,F) 
in law, in C(M+) x C(M) x C(M). 

Proof:. We first choose two arbitrary functionals <F e C b (C(R) xC(R)) and 6> G C b (C(R+)). We have 
&{W £ {a),F £ {a))G{X e > x {u,o-))P(da)W'' T (dw) 



P(da) 



V(W £ (a),F e (a))[ / 0(X £ ' x (u;,a))P^(du;) )P(da) 



*(W £ (a),F £ (a)) 



(0(X £ - x {uj, a)) - e{X x {u)))W' a {duj) 



Theorem 2.1 tells us that 



Q{X X {uj))f e - a {duj) x / <F(W e (a),F £ (o-))P(d<j). 



6»(X £ ' x (w, C r))P £ ' [T (dcj)^ / e(X x (w))F(duj), 

P a.s., and from Lemma 3.6, 

¥(W e {o),F e (v))P(dtr)-> f *(W{v), F(a))P(da), 
: J S 

as e — > 0. Hence, from the Bounded Convergence Theorem, we conclude that 
>F(W e (a),F e {a))e(X £ ' x (uj,a))P(d(j)P E '' y {doj) 

C 0{X X {u))P{duj) X / «P(W(<7),F(<7))P(d<7). 



It now suffices to note that A := {\P <& 0, & € C&(C(K) x C(K)), e C&(C(M+))} is a determining class on 

C(K) x C(M) x C(K+). □ 

It follows from Lemma 3.8 that 

^P) = W!) +; 1 /V^^ds + M/'*, (3.6) 



e y \ £ / 2-y/e 

in other words 

F e (xr)=F e w+^ t £ ' x +^r J 
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where {M^' ,t > 0} is the continuous martingale 



M^ x 



Jo 



X E ' X 



dM, 



and M x ' denotes again the martingale part of the process X E - X . In particular the quadratic variation of 
M e ' x is given by the quantity 



(M E ' X 



e I a 



X E -' 



1 



o a{Xf x /e) 



X e ' x 



£ 

Xf' x 



ds 



XI 



n 2 



ds 



\r—AW E (X^)) 2 ds, 



/o a{Xf x /e) 

and the joint quadratic variation of M E ' X and Z E is 



(M e ' X ,Z e ) t = y/k~ / & 



X E ' X 



i+x' 



X?<< 



V 



ds 



o a{XP x le) 



ds. 



Lemma 3.8. The identity (3.6) holds Q e a.s. 

Proof. Let {X%' X ' M , t>0} denote the process {X^' x , t> 0}, killed when exiting the interval [— e(M + 
l),e(M + 1)], and L £ ' a ' M its generator. For any M > 0, the random function 



<2>mO) 



Jo" ^k) So CM (f ) d ^ dz ' if ^ > °; 
/° /° c « (f) d y dz ' if x < °> 



with the random field {cm(x), —M — 1 < x < M + 1} defined by 

r /3m , if -Af - 1 < x < -M ; 
c M (x) = lc(x), if -M<x<M; 

Va-M, i£M<x<M + l, 



where 

a M 

Pm ■ 
satisfies 



I M+1 (i/ajz)) Iq AM c (v) dy dz ^ 
J™ + \z-M)/a(z)dz 

f-M-A l l a ( z )) fzv-M c (y) d y dz 

j:™_ 1 (z + M)/a(z)dz 



'(L E ^M)(f)=£- 2 CM(|) ! 

£m(-M-1)=<?m(M + 1)=0. 
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Consequently for each e > 0, 



and by virtue of the Ito-Fukushima decomposition (see [7] or Theorem 0.10 in [11]), we get for \x\ < M + 1 



*A-rf 



where 

r|f =inf{«>0,|Xri=<M + l)}. 
The result follows, by letting M — > oo, with the help of Lemma 3.2. □ 
Define 

h e (t): 



lo a(XP x /e) 
Lemma 3.9. As e — ► 0, 



a 



m C(R+) in probability. 

Proof. Denote 9{x) = — A. Then 9(x) is a bounded stationary field with zero mean. Letting 

e{x) =[{iF)l? [y)iv ) 

and repeating the argument in Lemma 3.8, we get 



t / ve.x \ / rXg' x /e \ 2 

-1 / 



where 

In the same way as in the proof of Proposition 3.7, one can show that the families {e 3//2 (9( — — ) — e 3 / 2 (9(-)} 
and {.Mf} are tight in C(R). Indeed, is constructed from exactly as <P from c. Moreover, 8 is, exactly 
as c, a stationary mixing bounded and zero mean random field. Now, multiplying the relation (3.7) by y/e, 
we conclude that 

ye.x \ 



in C(R+) in probability. This implies the desired convergence. □ 
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For any / G Cfe(R) and e > 0, let us define the process {N/' £ ;t> 0} by 
■dMf 



lo a{XT/e) 

rt f{xr)dzi 

We now prove the following proposition. 
Proposition 3.10. P a.s., 

{N } ' e ,X £ > x )=> (N f ,X x ) m C([0,i\) x C([0,i]), 

where 

Nf = l f f(X x )dX x . 
a Jo 

Proof. Since (X e - X — Z e ) converges to zero in probability (N/' e , X^' x ) behaves as e — * exactly as 
(iV/' E , ), hence we consider the two-dimensional martingale (iV/' E , Zf), and compute its associated bracket 
process, which takes values in the set of 2 x 2 symmetric matrices. We have 

(Nt>*\\\ _( So a'i?/s) dS ~ a ti c!{xf/e) ds \ 
\ ///* \ a Jo a(^'*/e) dS a JO a^'VO / 

Combining Theorem 2.1, Lemmas 3.5 and 3.9 we obtain that this R 4 -valucd process converges P a.s. in P 
law towards 

ljlf{X x )ds J*f(Xf)da 
ftfiXfida at 

We then conclude that P a.s., 

(N f ' e ',X B ' x )=>(N f ,X x ) 

in P law, where 
i- 1 



N( = l f f(Xf)dXf, t>0. 
a Jo 



□ 



The statement below is a straightforward consequence of Birkhoff 's crgodic theorem. 
Proposition 3.11. For any N > 0, iel and f € C(R) the following convergence holds P a.s. 



sup 

\y-x\<N 



Jx a{z/£) aj x 



0. 



We now establish the version of (3.6) for e = 0, which is an Ito type formula for the process {F(x + Xt),t > 
0}, where F(y) := £ J Q y W(z) dz, ygl, More precisely, we have the following lemma. 

Lemma 3.12. For any t>0 and ifl, we have 

F(X X ) - F(x) = I f W(X X ) dX s + 1 / Lr x W(dy). (3.8) 
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Proof. Wc prove this formula by using smooth approximations of the process {W^}, obtained by convolution. 
Let p be a Cq (R) function such that p > 0, supp(p) C [—1,1] and L p(y) dy = 1. Define now 



p„(y):=np(ny), W n (y) := p n (y - z)W(z) dz = p(z)W[y--)dz 



(3.9) 



From the uniform continuity of W on compacts, ||Wn — W||c(jf) —* 0, P a.s., as n — ► oo, for any compact set 
A' in R. Moreover, taking into account the fact that W is a standard Brownian motion, we get 



E(W*(y)) < C e(w 4 (y-ffj ^ = C £ (y - £) dz < C(l + y 2 ), 



(3.10) 



for y e R. Set 



F n (-):=- / VF„(y)dy. 
Ito's formula applied to the process {F n (x + X t ), t > 0} gives 



F n {x + X t ) - F n (x) 



W n (x + X s )dX s + - I W' n {x + X s )ds. 



(3.11) 



Recall that {x + X t , t > 0} is a non-standard Brownian motion independent of {W}. It is easy to see that 
the left-hand side in the last formula tends to F(x + X t ) — F(x), P x P a.s., as n — > oo. Moreover, 



E x E 



(W n (x + X s )-M/(a; + X s ))dX 
(W„(x + X s )-iy(x + X s )) 2 ds 
< afij^ [E(W n (x + X s ) - W(x + X s )f] ds 



= aE{ E 



->o, 

since W n (x + X S ) — > W(x + X S ), ds x P x P a.e., as n — > oo, and moreover the sequence {(M /r „(a; + ^ s )) 2 , n > 
1} is ds x P x P uniformly integrable on [0,t] x x J?, thanks to (3.10). 

Finally, from the occupation time formula for continuous semimartingales (see e.g. Corollary 1.6, page 
209 in [15]), with {L' t } denoting the the local time of the process {X s , < s < t}, 



fw' n {x + X s )ds= I Lr x W^(y)dy^ f Lt x W{dy), 

JO JR JR 



in L 2 (S), P a.s., as n — > oo (for more details see Section 5.7 in [10]). We used again the fact that the 
Brownian motions {X t ,t > 0} and {W(y), y S R} are independent. Passing now to the limit in the formula 
(3.11) we get the desired result. □ 

We can finally proceed with the following proof. 

Proof of Theorem 3.1. Since the mapping 

A:C(R+) x C(R)^C(R+) 
defined by 

A(x,f)(t) = f(x(t)) 



Homogenization of random PDE 



5.31 



is continuous, if we equip the three spaces with the topology of uniform convergence on compact sets, we 
first conclude from Proposition 3.7 that 

W £ {Xf' x )=>W(X x ). 

Hence from the formulas for (M £ ' x ) t and (M e ' x , Z £ ) t above, and Lemma 3.9, we deduce easily that 



(M e ' x ). =>> — / W(X x )ds, 
a Jo 

(M £,X ,Z £ ), => c f W(X x )ds, 
Jo 

and consequently 
a Jo 

From Proposition 3.7, those convergences are joint with those of (Xf' x ,F e ). Consequently 

F e (Xn - F e (x) ~ Mt x => F(X X ) F(x) -if W(X X ) dX s . 

a Jo 

The convergence Y t e ' x => Y x now follows from (3.6) and (3.8). The result finally follows from the fact that 
all the above convergences are joint with that of Xf' x . □ 

4. Pointwise convergence of the sequence u E 

The first part of this section is devoted to establishing uniform intcgrability estimates for the exponent 
in the Feynman-Kac formula (Propositions 4.4 and 4.5) which are essential for the proof of the pointwise 
convergence part of Theorem 2.2, to which the second part of this section is devoted. 
We first define the following R+ -valued random variables, for 0<7<l/2, e>0: 

c \W s (x)\ 
£ 7 , e = sup ■ 



We have the following lemma. 

Lemma 4.1. For any < 7 < 1/2 and So > 0, the collection of random variables {£ 7 , e , < e < £0} is tight. 
Proof. Due to the symmetry it is sufficient to estimate |W e (x)| for x > 0. We have 

E(\W £ {r)\ 2 ) = s f ^ f ^ E(c( S )c(t))dsdt 
Jo Jo 

rr/e 1*00 

<2e / \E(c(0)c(s))\dsdt 



Jo Jo 
< 2rc a . 

Denote 

/•OO 

r) t = / E(c(s + t)\g t )ds. 
Jo 
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By Proposition 7.2.6. in [6] the process rj t is stationary and \rj t \ < c\ a.s. with a non-random constant C\. 
Moreover, 

t 

c(r) dr — r/t 



is a square integrable £/t martingale. Denote it by TV*. Clearly, 
W e (f) = ^ /' ' c(s) ds = ^N t/e + ^ %/£) 

C Jo c c 

and thus we deduce from Doob's inequality 



E( sup |V^(i)| 2 ) < ±E( sup (V~eM t f) +2^f 

M)<t<r 7 C \Kt<r/e 7 C 



<4^((^ /£ ) 2 ) + 2$ 
cr c 



<8E(|W e (r)| a ) + 10^ 
<C(e + r), 

provided C = (16c ) V (10cf/c 2 ). Now for j >1,M> 0, 
\W E (r)\ 



P( sup J^V' >M )<P( sup |^ £ (r)|>(l + 2^ 1 ) 1 ^M) 



< 



0<r<23 

C{e + V) 



M 2( 1 + 2 J-l)2-2 7 

<( £ Vl)|^(l + 2^) 2 ''- 1 . 



Summing up over j > 1, we deduce that 

P(e 7)£ >M)<2p(sup I i^ ? >M 

^-Vl)^E(l + 2^- 1 



j'=o 

A/ 2 



The lemma is established. □ 

Remark ^.2. We can in fact show that, as e — > 0, 

\W{x)\ 
xeH (1 + Fl) 1 t 

provided again < 7 < 1/2, 6m< we shall not use that result. 

We next state a result, which is an immediate consequence of Lemma 3.2. 
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Lemma 4.3. There exists a continuous mapping 

p-M+ xR+ x (0,2)^R+ 
such that for all c> 0, t > 0, e>0, 0<p<2, 

W'-cxp(c\X^ x \ p )< p ( c ,t,p). 

We next establish the following proposition. 
Proposition 4.4. For any < 7 < 1/2, there exists a continuous mapping Wj :K + 
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such that 



cxp 



1 



ds 



(4.1) 



Proof. Since from (3.6) and (3.5) 



1 



X e ' x 



we obtain 



ds = 2 



a{y/e) 



Ay - Ml 



W'' exp 



1 



< E e ''exp 8c 



X e ' x 



ds) =E e '-exp 4c 



1/2 



V W s {y) 
a(y/e) 



dy - AMI 



W £ (y) 



dy (E e '-exp(-8Aff)) 



1/2 



lx a(y/e) 

Clearly, it suffices to estimate each factor on the r.h.s. of (4.2) separately. 



(4.2) 



E e -- exp ^8c 



-V; 



" W e (y) 
a(y/e) 



dy < E e -' exp 



< E e '' exp 



A'f 



x; 



\W s (y)\dy 



^(i + lyD^dy 



< E - exp (c'^[(i + \xr\f-~< - (1 + M) 2 " 7 ; 

where we have used Lemma 3.2 for the last inequality. The second factor on the r.h.s. of (4.2) can be 
estimated as follows 



E e -' exp(-8Mf £ ) < (E £ -- exp(-16M t £ - 128(A/ e ) f )) 1/2 (E e '' exp(128(M £ ) t )) J 

The first term on the r.h.s. does not exceed 1. For the second one we have by Jensen's inequality 

ft 



,1/2 



E £ '-cxp(128(A/ £ ) t ) <E £ ''exp^c"^ W*(X?*)da 

<i -1 / E E '-cxp(c"tWZ(X £ s > x ))ds 
Jo 



< sup E^exp(c^K 7 , £ ] 2 (l + |X|^|) 2 - 27 ) 



0<s<t 
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where we have again used Lemma 3.2 for the last inequality. The result clearly follows. 
Clearly, the same proof allows us to establish the slightly more general proposition. 



□ 



Proposition 4.5. Let {t £ , e > 0} be a collection of stopping times such that < r e <t PxP a.s. Then for 
any < 7 < 1/2 , 



exp -= 



1 rv*£ 

V £ 



ds 



where 



R + is the mapping which appeared in (4.1). 



We can now proceed with the following proof. 

Proof of the pointwise convergence in Theorem 2.2. We will now show that for each (t,x) £ R+ x R, 
u £ (t,x) => u(t,x), as e — •> 0. We delete the parameters t and x. It suffices to show that for any tp £ C(R; [0, 1]), 
Lipschitz continuous, as e — * 0, 



where X £ = X e t ' x , X = Xf = x + Va~B t , and (recall (3.6) and (3.5)) 



(4.3) 



Y = 2- 



X' 



w £ {y) 

a{y/e) 



dy - Ml 



2c 



X; 



" WM dy 1 

a(y/e) a 



W £ (X e s ' x )dZl 



x 



W{y)dy- / W(x + X s )dX s 



The fact that this Y equals the exponent in the Feynman-Kac formula for u(t,x) follows from (3.8). We 
first approximate Y e by Y e - M as follows. For each e > 0, M > 0, let 



r e M =M{s>0;\XI' x \>M} 



and 



a(y/e) a 



W e (X e s x )dZ\ 



We postpone the proof of the following lemma. 
Lemma 4.6. 

sup|£V(E e <- [g(Xn cxp(F E )]) - E(p(E £ '' [g(X s t < x ) exp(F E ' M )])| - 0, 

e>0 

as M — ■+ 00 . 

Since the collection of random processes {W e (y);y £ R} is P-tight, for all S > 0, there exists N £ N and 
/<5,i, /<y,2j ■ • • , /«,jv e C 6 (E + ) such that if 



:={ sup \W e (x)-f 8 ,k(x)\<s}, 



Kk< N, 



then 



p((f>)>* 
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Let now 



and for 2 < k < N, 



fc-i 



1=1 

and finally 



so that P(A S ' e ) < <5. Note that ||/<5,fc||oo depends on 6. However, we can and will assume that for some 
0< 7 <!/2, 



|/ 5 , fe 0r)|<2£ 7 (l + |z|) 

for all 8 > 0, k e N. 
We now develop 



1-7 



(4.4) 



JV 



i^(E £ >-[ 5 (X e )e yt ' ]) = ^£;{^-[ 3 (X e )e ye ' ]);^ £ } + J B{^(E £ '-[ 5 (X e )e yt ' 

k=l 

The last term in the above right-hand side is bounded in absolute value by 6. Now for 1 < k < N, 



£{^(E S '" [g{X^° ']); B s k > s } = E{<p(E £ '' \g(X e )e Y **]);B%.' e } + ef , 



where 



Y°>* := 2c 



/ 4 , fc (Xf a =)dZ] 



a(y/e) 5 J Q 
We postpone the proofs of following lemmas. 

Lemma 4.7. There exists a constant C , which depends only on t, \\fs,k\\oo and the constants appearing in 
(2.1), such that for each 6 > 0, 1 < k < N, a G E, 

supE £ '' T [e 2 ^' e ] <C, P a.s. 
Lemma 4.8. 

N 



sup 

e>0 



E 1 

k=l 



0, 



as 5-^0. 

It follows readily from Propositions 3.10 and 3.11 together with Lemma 4.7 that, as e — > 0, 



E e -<W) e V ]^E^[ ff (X)e^], 
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P a.s., where 



A';' 



fsAv)*y- / f6,k(xf)dx 3 s 



Let E>1 , 1 < k < N, denote the sets defined exactly as the B k ' £ 's, but with W e replaced by W. The boundaries 
of those sets being of zero Wiener measure, we conclude from the last statement and the fact that W e =>■ W 
that as e — ► 0, 



i?{^(r-[ g (r) C ^]);B^}^i?{^-[ S (I)e^]);B^}. 
Now, in the same way as above we obtain 

N 

^(E^-[ 5 (X)c^) = ^S{^-Lg(X) e y ]);^} + S{^-L 9 (X)e^]);^n, 
fe=i 

P{A S ) < 6, and for each 1 < k < N, 

E{ V (W>- [g(X)c Y })- Bi} = E{y{W>- [g(X)e^]);B s k } 
All we need to conclude the proof is the next lemma. 



el 



Lemma 4.9 

JV 



E 



k=l 

as 5 — ^ . 

It remains to prove the four lemmas. 



Proof of Lemma 4.6. Let K denote the product of the Lipschitz constant of ip by |j.g||oo- Since (p takes 
values in [0,1], 

E\p(E^[g(X s )e YC }) <p(E^ [giX^e^])] 
<E{lAKE^(\c Y ' -e yE ' M |)} 

< E{1 A V2K(W'-[ C 2Y ' + e 2Ye ' M )) 1/2 (F s >-(r e M < t)f 2 } 

< E{1 A [2K (f 7 (^, e )) 1/2 (P^(rI / < i)) 1/2 ]} 

< L(^-(T S M < i)) 1/2 + P(*y(^,e) > L 2 /4K 2 ), 

where < 7 < 1/2 and !^ 7 (£ 7ie ) appears in Propositions 4.4 and 4.5. Since the latter is P-tight (for e > 0), 
we can choose L = L n such that the second term in the last expression above is less than rj/2. It remains to 
choose M large enough (exploiting this time the tightness of {X £ ' x ;e > 0}) such that 



^(P £ --(r| f <t)) 1/2 <V2. 



□ 



Proof of Lemma 4.7. There exist two constants C\ and c 2 such that 
E^c 2 ^" < (E e --c c ^*"-l) 1/2 (E^c C2 /o fsMKn^y/\ 
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The bound for the first factor on the right follows from Lemma 3.2, and the bound for the second factor 



follows easily from the boundedness of both fs y k and 4-(Z £ ^ 



Proof of Lemma 4.8. We have, by an argument similar to that in the proof of Lemma 4.6, 



N 



N 



k=l 



5>£ 4 \<E lA^KE^e* 

V k=l 
( N 



■e^)|y^_^ )£|]li 



Now 



fc=i 



< (2E--[c 2ye ' 5 " +c 2Y «"]) 1/2 (l K ..W>-{\Y^ 1 -Y^\ 2 }) 1/2 
<c(2[f 7 (£ 7 , £ ) + C]) 1/2 



X(W(I«« \ 9 )+^'(\xr-X^ \ 2 )+5t + W'(t-tATt-i)) 1/2 l B », 



Finally. 



iV 



£|ef | < £(1 A [#(1 + ^, e )f\5t + CS 2 + p(t, <5)) 1/2 ]), 



k=l 



where 



p(t,S) =sup[E £ <-(LYr-^ j'j+E^Ct-tATf-O] 



£>0 

0. 



5-1 



as S — > 0. The end of the proof is similar to that of Lemma 4.6. 
Proof of Lemma 4.9. This proof is similar to that of Lemma 4.8. 



/v 



JV 



Y,E{\^ e -[9(X)e Y }) - V (W-lg(X)c Y » })\; B s k } < e\ 1 A K^E e '-\e Y - e Y * \l B s 



k=l 



fe=l 



Now 



l Bl E--\e Y - | < (2W>-(e 2Y + e 2 ^)) 1/2 (l fl jE e - [|Y - Y k s f}) 1/2 , 

k k 

and on the set B k , using in particular (4.4), 



E e '-(i^-n*i 2 )<8( ~ 



[^(y)-Afc(y)]dy 



+8 UJ E£ ' 7o i^+^-^^+^r^ 

< cUH + ^t 1 ' 2 V ^-^^.-(iBtl > S- 1 -\x\) 
+ c£ 2 t 2 -<F £ '-( sup Ix + BsIxT 1 " 



□ 



□ 
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which clearly goes to 0, as 8 — > 0. The result follows. □ 

□ 

5. Convergence in C(R_|_ X R) 

It remains both to prove the convergence of the finite dimensional distributions of u E towards those of u, 
and to establish that the sequence {u £ ; e > 0} is tight as a collection of random elements of C(R+ x R). 

Theorem 5.1. For any £ e N, (ti,Xi) el + x R, l<i<£, 

{u E (ti, x{), u e (t t , xe.)) =>• (u(ti,Xi), . . . , u(U, x e )) 
as e — > . 

Proof. We only sketch the proof, the details being identical to those of the proof of the pointwise conver- 
gence, as were given in the previous section. For each 1 < i < £, wc define Xf := X^ Xt and 



o \ e 

We need to take the limit as e — -> in the quantity 

^(E £ >-[ 5 (Xf)exp(Ff)], . . . ,E s <-[g(Xt) exp(Y/)]), 

where ip <G C(R f ; [0, 1]) is Lipschitz continuous. For that sake, referring to the notations in the previous 
section, for each 6 > 0, 1 < k < N , wc define 



Xf i r\ 1 rt 



where 

/ Y e ' : 

z t =X t + ex I — 
We have that for each 1 < i < £, 

N 

[g(X!) exp(Ff )] ~ £ E e '' [ 5 (Xf ) exp(F^)] 1 B , , 
fc=i 

and consequently 

B V (E e '-[ff(Jff) exp(F 1 e )], . . . ,E e --[ 3 (Xf ) exp(y/)]) 

N 

= J2E[^[g(X!)exp(Ylt)l ■ ■ . ,E-[. 9 (X|)exp(r/,:)])l B ,]. 
fc=l 

But for each 1 < i < £, 1 < fc < N, 

[g(Xf) cxp(y^)] -> E £ >" \g{Xi) exp(y^)] 

— > 0, where X; = x, + \/hBt i and 
lif fc = 2 1 (J * dz - J " /j ifc (xi + v 7 ^) dB f 
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The result follows. 

We finally establish the tightness result. All we need to show is the following theorem. 

Theorem 5.2. There exists a sequence p £ in R + which tends to as s — > 0, a mapping 

(f : R\ x R 2 -> R + 
such that for all T > 0, M > 0, 

sup (p(T,M,s,t,x,y)^>Q 

0<s,t<T;\x\,\y\<M;\t-s\<S;\x-y\<8 

as 5 — > 0, ond a continuous mapping <P : — > R + swc/i i/iai /or all < s,t <T, \x\,\y \ < M , 

\u(t, x) ~u(s,y)\< <P(T, M, f 7>e ) [p e + V (T, M,s,t,x, y)]. 
Proof. We will establish (5.2) with 

¥>(T, M, s, t, x, y) = y/\t - s\ + y/\x -y\ + ^E([g(x + Va~B t ) - g(y + Va~B s )] 2 ), 

which satisfies (5.1). 
We need to consider 

u'(t,x) - u s (s,y) = E £ <- \g(X e t ' x )e 1/VE fo< x °'*/ £)dr ] -W \g(X^)e 1/VI f» ciK ' y /e)dr 



'■{[g(X £ t ' x ) - g{Xl>y)]e 1,V ~ e $o< x ^/^ A ^ 



l/VZ J* c(X^/e) dr _ 1/V5 J° cW/e) dr 



The absolute value of the first term of the right-hand side of (5.3) is dominated by 



\jw^ 2, ^^ c{X ^ /e)dr x ^W<-{\g(Xt> x )-g{XV v )] 2 ). 
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(5.1) 



(5.2) 



(5.3) 



The first factor in the last expression contributes to the coefficient <?(T, M, £ 7 , e ) in (5.2), while the difference 
between the second factor and 



y/E e <-(\9(x + V&Bt) - g(y + Va~B s )} 2 ) 



contributes to p E . Now the absolute value of the second term in the right-hand side of (5.3) is dominated 
by cE(A e +B £ ), where 



and 



.4, 



Br 



l/y/S f* c(X^/e) dr _ 1/y/Z f* c(Xpy/e) dr 



l/Ve f o c{X' r -y/e)dr J* c(X°' * /e) dr 



Below we will use repeatedly the elementary inequality \e a — e b \ < \a — b\e aVb . In particular, we deduce that 
(assuming w.l.o.g. that < s < t) 



Bp < 



Xi'V 



dr 
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Consequently EB £ is dominated by a factor which contributes to #(T, M, £ 7jE ) times the square root of 



27 



X; 



W £ (z) 
a{z/e) 



-I 2 



dz + Mf' v - M £ > y 



<<Jt-s\ 



sup W 2 (z)\X £ t ' v -X%*\ 2 + [ W 2 {X*' y )dA 
'• y \y\xl- y \ Js J 

(i + \x e r >y\f- 2 ~<dr 



i+ sup \xpy\f ^\Xf v -X^f 



s<r<t 



By the same argument as above, W''A £ is dominated by a factor which contributes to #(T, M, £ 7 . e ) times 
the square root of 



x 



'" W e (z) 
a{z/e) 



dz 



xf<» 



W £ (z) 
a{z/e) 



dz 



\Ml>*-Ml 



(5.4) 



While the two laws of {Xf' x } and {Xf' y } are given to us, the coupling between these two processes is at our 
disposal. We make the following choice. Given a € S, the two processes {Xf' x } and {Xf ' v } arc mutually 
independent, until the first time r^ y when they meet, and then the two processes follow the same trajectory. 
We note that as e — > 0, 

T xy =^ T xy = hif{r; x + VhE>l = y + Vfi-B 2 }, 

where {B 1 } and {B 2 } are two mutually independent standard Brownian motions. This follows from the fact 
that the first time when a two-dimensional Brownian motion meets the diagonal of R 2 is a.s. a continuous 
function of the trajectory. Suppose w.l.o.g. that y < x. Then 



■\T xy >b) 



= 1-: 



sup (B 2 r -B l r )<^l 
o<r<5 ya 



>•( sup (B>-fli)>f» 
\0<r<s ya 



= 1 - 2P £ '- 1 Bj - Bg > 



i^ x -y 



z\< 



1 Ire — 7/1 

where Z is an N(0, 1) r.v. On the other hand, 

E £ '' (r xy At) < c\fi\x — y\. 
We now estimate the first term in (5.4) for |a;|, \y\ < M, < t < T, 



x 



< 



°" W e (z) 
a{z/e) 

W e (z) . 
a(zje) 



dz 



x? 



dz 



W £ (z) 
a(z/e) 

Xr W e {z) 
a(z/s) 



dz + 



X; 



W £ (z) 
a[z/e) 



1 {r£y>t} 
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< ^M)^]* - y| 2 + E-v[(i + | Z | v v \y\ v |^|)^i {xSy>t} ]) 

<c(M,r)£y|z-y| 2 + |z-y|+py, 
where p' e — > 0, as £ — > 0. We finally estimate the second term in (5.4), again for \x\, \y\ < M, < t < T, 



E^(\Mt x - Mt y \ 2 ) < cW r yM [W?(X^) + W?(Xfv)]ds 

Jo 



sup (i + \xr\ + \xt'y\f- 2 \ T i y M) 



0<s<t 



<c(M,T)e,J\x-y\+ P ' £ ). 



The theorem is established. 



□ 



6. The stochastic PDE for the limit u 

In this section, we study the limiting SPDE. Formally, we would expect that it reads 

^(t,i)=2^2 (*> x ) + 2u (*> x ) W(dx), t>0,xeR; 
u(0,x) = g(x), xGM., 



(6.1) 



where the above stochastic integral should be understood in the sense of the anticipative Stratonovich 
integral (see [13] or [12]). However, since it is difficult to get any uniqueness result for such an equation, we 
prefer to rewrite it in a different form. In fact, since 



1 fx 

Je \e 



_dW £ 



we can rewrite the original u e equation as 

%^ X) = \§x SO +^(^<0(M) -W. 

u e (0,x) =g(x). 
Hence we might expect that the limiting equation reads 



^(t,x) = -^{t,x)+c—(Wu){t,x)-cW(x) — (t,x), t>0,x&. 
u(0,x) = g(x), i£i 



(6.2) 



Would the equation be posed on a compact interval, or would the random process {W(x),x G K} have a.s. 
bounded trajectories, then the existence and the uniqueness for (6.2) would be very easy to obtain. Here we 
will content ourselves with the following result. 

Theorem 6.1. The parabolic PDE (6.2) has a solution ue Lj 2 oc (R+; H^ oc (M)) a.s., which is given by the 
Feynman-Kac formula 



u(t, x) := . 



ff (Xf)«p(| / Lr x W(dy) 



(6.3) 
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Proof. We need to show that u, given by the formula 



u(t, x) = E 



fl(Xf)«p(| / L*- x W(dy) 



belongs to Lf oc (R + ; Hl oc (R) ) , and solves the parabolic PDE (6.2). 

For that sake, we first define an approximation of the Wiener process W, by the formula 

W n (x) = [(W*p n )(x) An] V(-n), 

where * stands for the convolution operation, and p n {x) := np(nx), p is a smooth map from R into R+ with 
compact support, such that 



/ p{x) Ax = 1. 

JUL 



Let u n denote the solution of the approximating PDE 



(t,x) + cW„(x)u n (t,x), (>0,i£M: 



dt ^ X> ~ 2 cb 2 
u n (0,x)=g(x), xeR. 

It follows from the Feynman-Kac formula, see e.g. [9], that 



(6.4) 



u n (t,x) =E 



E 



, 9 (.x + JT t ) exp c / W'Jx + X s ) ds 



g(x + X t )exp( - / W^z)L z - x dz 



(6.5) 



where X t = V&B t , and {5} is a standard Brownian motion defined on (Q, T , P), while Lf( denotes its local 
time at time s and point y. It follows from arguments similar to but simpler than those in Section 4 that 



E 



g(x + X t )exp[^ / W' n {z)Ll- x dz 



E 



g{x + X t )ey^ \ - / Lf-*W(dz) 



For each M > 0, we now write an equation satisfied by u r 



~ad 2 



du 11 u 



(t,x) + c 



(t,x)-cW n (x)^—(t,x), t>0, -M<x<M 
ox ox 



u n (0,x)=g(x), u n (t,-M) = C(t,-M), u n (t, M) = £ n (t, M), 
where £" denotes the right-hand side of (6.5). It is now easy to show that 

£ n (t, M) - £ n (t, -M) £"(*, M) + f (i, -M) 



(6.6) 



v n (t,x) :=u n (t,x)-x- 



2M 



solves eq. (6.6) but with homogeneous Dirichlet boundary conditions. Now v n converges strongly in L 2 oc (R + ; 
Hq(—M,M)), P a.s., towards the solution of the parabolic PDE 

^ (t,x) = 1 (t, x) + (t,x)-cW(x)^(t,x), t > 0, -M < x < M; 

v(0,x)=g(x), -M<x<M, v(t, —M) = v(t, M) = 0. 
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We conclude the u := linin^oo u n belongs to the space Lf oc (R+;iJj] 3C (K)) a.s., and it satisfies (6.2) in the 

variational sense, i.e. for any t > 0, any ip e C 2 (R) with compact support, and a.s. ((•,•) denotes the scalar 

■2/ 



product in L 2 (R)), 



(u(t),ip) = (g,cp) + 



-{u(s),cp")-c(Wu(s),ip')-c 



□ 



Remark 6.2. The same problem on a bounded interval (a, b) Cl with Dirichlet homogeneous boundary 
conditions could be treated either by the same method, or by a much simpler PDE argument. Namely, noting 
that 77gc(^) = cW' e {x), we have that u e is the unique element of Lf oc (M. + ; Hg(a, b)) H C(R+; i 2 (a, b)) which 
satisfies the PDE 

(¥ = ^(«(DIf) + ^(^ £ )-^^, a<x<b, t>0; 
\ u £ (0, x) = g(x), u e (t, a) = u e (t, b) = 0. 

It is then not very difficult to show that u £ =>■ u, the unique solution again in L 2 oc (R + ; HQ(a,b)) (~1 
C(M+; L 2 (a, b)) of eq. (6.2) in R + x (a, b), with homogeneous Dirichlet boundary conditions. 
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